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Abstract 

We present a theory of the melting of elemental solids as a dislocation-mediated 
phase transition. We model dislocations near melt as non-interacting closed strings 
on a lattice. In this framework we derive simple expressions for the melting tem- 
perature and latent heat of fusion that depend on the dislocation density at melt. 
We use experimental data for more than half the elements in the Periodic Table to 
determine the dislocation density from both relations. Melting temperatures yield 
a dislocation density of (0.61 it 0.20) in good agreement with the density ob- 
tained from latent heats, (0.66 it 0.11) b^^, where b is the length of the smallest 
perfect-dislocation Burgers vector. Melting corresponds to the situation where, on 
average, half of the atoms are within a dislocation core. 
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1 Introduction 

Nearly 50 years ago Shockley JIl successfully accounted for the fluidity of a liquid by 
assuming a certain concentration of line defects in the liquid state. Bragg [0 had earlier 
estimated an upper bound on the core energy of a dislocation under the assumption 
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that the core atomic configuration was hke that of a hquid. Cotterill and Doyama 
later confirmed the Bragg estimate. These early results implied that the liquid state is 
equivalent to a crystal saturated with dislocation cores. It was first suggested by Mott 0] 
that the melting transition could be described in terms of dislocations. 

There are now compelling results from molecular dynamics and Monte Carlo 
simulations that imply that dislocations play a key role in three-dimensional melting, and 
moreover there is experimental evidence that linear defects are in fact generated near the 
melting transition |]^. 

Mizushima ^ and Ookawa ^ were the first to formulate a dislocation theory of melt- 
ing. They based their theory on the fact that the self-energy of a dislocation decreases 
with dislocation density because of screening. In their theory melting is a first-order 
transition that occurs when the free energy of the crystal with a sufficiently high concen- 
tration of thermally-generated dislocations equals the free energy of the dislocation-free 
crystal. Their predicted melting temperatures agree with data for reasonable choices of 
the core energy. Many dislocation theories of melting then followed [10-13], most notably 
the exhaustive treatment of linear-defect-mediated melting by Kleinert and collaborators 
ril. We refer the reader to several fine reviews of the literature for additional details and 



references on dislocation-mediated melting |15, 16 



Significant progress in our understanding of melting has been achieved by Kleinert 
iM who pointed out that the melting process cannot proceed through the mediation of 



dislocations alone. Dislocations are associated with the discrete translational symmetry of 
the crystal, so only this symmetry is lost when dislocations condense. But the rotational 
order of the solid is also lost as the solid converts into liquid, and for this to occur the 
defects associated with the rotational symmetry of the lattice, namely disclinations, must 
come into play. Kleinert assumes that the free energy of dislocations alone would lead to 
a second-order phase transition, and in addition shows that a second-order proliferation 
of disclinations can occur in a background of dislocations above some critical density. The 
coupled dislocation-disclination system could, however, undergo a first-order transition, 
i.e., melting. 

Although disclinations must participate in the melting process, in this paper we con- 
sider only the dislocation degrees of freedom. Ideally, we would derive the precise form of 
the free energy of a dense ensemble of dislocations interacting via the full Blin potential 
1^ and subject to specific configurational constraints (Brownian, self-avoiding, etc.), but 



this problem has so far defied solution. Instead we develop an effective theory of melting 
based on perfectly screened, non-interacting dislocations. We employ the widely accepted 
—p In p form (p is dislocation density) for the self-energy density of dislocations p, |1^, |T^ , 
which results in a first-order phase transition. A dislocation in a dense ensemble of other 
dislocations is assumed to be a random loop, i.e., the possible configurations of a disloca- 
tion loop are closed random walks, and short-range steric interactions are neglected. Thus 
the partition function is evaluated in the independent-loop approximation. We obtain two 
new relations: a simple expression for the melting temperature (the melting relation) that 
explicitly takes into account the crystal structure, and another relation between melting 
temperature, latent heat of fusion, and critical density of dislocations. We carry out a 
comprehensive comparison of these relations with experimental data on over half of the 
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Periodic Table. The melting relation is accurate to 17% [|r^. Dislocation densities as de- 
termined from the melting temperature and latent heat relations are p = (0.61 ± 0.20) 
and (0.66 ±0.11)6"^, respectively, where b is the length of the smallest perfect-dislocation 
Burgers vector. Both relations should also apply to alloys and compounds. 

In Section 2 we discuss the statistical mechanics of dislocation loops on a lattice. 
These results are used to derive the melting relation in Section 3, the free energy density 
in Section 4, and the formula for the latent heat of fusion in Section 5. The values for the 
critical dislocation density extracted from both the melting relation and the formula for 
the latent heat of fusion are checked in Section 6 with the formula for volume change at 
melt. Our concluding remarks appear in Section 7. 



2 Statistical mechanics of dislocation loops on a lat- 
tice 

The energy per unit length, a, of a dislocation can be very large. However, this energy can 
always be compensated at sufficiently high temperatures by the large entropy of line-like 
structures, as will be seen in what follows. 

In a Bravais lattice with coordination number z we consider the graph F, the edges of 
which are all nearest-neighbor links. The set of z links from any lattice site is identical to 
the set of shortest perfect-dislocation Burgers vectors, of length b. We now evaluate the 
partition function for a single Brownian loop on F. 

The line tension, a, is assumed to be independent of its length, L, as is the case 
for a dislocation in a dense complex network. (In a dilute network, interactions between 
distant segments of a dislocation lead to a logarithmic dependence of a on L.) The number 
of configurations of a string of length L is {z')^^'', where z' is the number of possible 
directions that a line segment can take from a given lattice site. If backtracking is not 
allowed, z' = z — 1. For a simple cubic lattice in D dimensions z = 2D. 

Hence, the partition function for a single closed dislocation (in 3 dimensions) is 

Z, = Y^piL, V){z'fl\-^^^ = Y^piL, V^)e-/^-«^ ^ a (^1 - ^) , (1) 

where (3 = l/fc^T, p{L, V) is the sum of probabilities over all lattice sites that a dislocation 
of length L will close, V is the volume of the system, and cXefr is the effective energy cost 
to create unit length of a string at temperature T. 

In order to calculate p{L, V), let p{r', r; L/b) be the probability density for a dislocation 
of L/b = n steps to start at r and end at r'. In the limit n oo, 6^0, and L = const. 



p(r', r; L/b) satisfies the diffusion equation |2y] 



,2 



= ^V^P. (2) 

on z 

the solution of which is the heat-kernel expansion 

p(r', r; L/b) ^ p{v' - r; L/b) = ^ ^(r)/*(r')e-^'<^/^ (3) 
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where the /k(r) are eigenfunctions of the Laplacian, 



z 

which we take as normahzed according to 



- -yh = i^k/k, (4) 



dh\Uv)\' = l. (5) 



1 /)2P 

h{r) = -y=e^^-, E^ = -^, (6) 



When — i> oo, we have 



where < = |k| < oo. It then follows from Eq. (3), in which we replace a sum by an 
integral, Ek ^ V/{2TT)^Jd% that 

Pir' - r; L/b) = / ^ ^Mr-r')-bLky.' ^ ^-.'(r'-r)V46L_ (7) 

The normalization (5) of the eigenfunctions thus imparts unit normalization to the prob- 
ability density: 

' d^up{u; L/b) 



The partition function for a dislocation loop (r' = r) is therefore 

(9) 

where the factor b/L removes the overcounting due to the degeneracy in the number of 
starting points on the loop. Here, N{L) is the number of configurations of a loop of length 
L. The exponent 5/2 becomes 1 + D/2 m D dimensions pT[| . 

Real dislocations are not necessarily Brownian loops. In fact, they are expected to 
be self- avoiding and/or neighbor-avoiding loops, so they do not penetrate each other's 
core. Eq. (9) can then be extended to non-Brownian or open dislocations by means of an 



effective exponent g + 1 7^ 5/2 and normalization constant A{q^z') [^, as follows: 



^i = ^(^>-O^E(f)""e-'"«" (10) 



Here, q = —1 for non-interacting (Brownian) open dislocations ||2^ and g ~ 7/4 for self- 
avoiding dislocations at low densities in 3 dimensions [^. In the string literature, the 
value q = has also been quoted. A general argument based on modular invariance 



shows that for non-interacting closed strings g = for sufficiently high energy on any 
compact target space. The same value of q was also obtained in a discrete model for 
strings [^, and as a static solution to the string Boltzmann equation ||2^. In principle, q 



may even be a function of temperature. Although we may expect 3/2<g<7/4 [|12|, our 
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main conclusions do not depend on the precise value of q. The normalization constant, 
A{q,z'), can be calculated analytically for Brownian loops in any dimension {q = D/2), 
analogously to the calculation of A{3/2,z') = {z'/AttY^'^ in Eq. (9), and numerically in 
other cases. 

The average length of a loop is 

_ Y.LLN{L)e-^'^^ _ e(T) 



where we define 



and 



e(r)^^(g,z')E(f) 'e-P^^-\ (12) 



L/h 

e(T) ^ A{q, z') Y: (f ) e-^<^-^ = ^ Z,. (13) 

L/b 

The grand canonical partition function for an ensemble of non-interacting indistin- 
guishable loops is given by 

Z = Z{T, V, /.) = E ^e'^^/'^-^ = exp{exp (^) Z,}, (14) 

where /j, is the chemical potential. The free energy of the ensemble is 

F = -keT In Z^ -ksTe^/'^'^Zi. (15) 
The average number of loops in the ensemble is 

N^-(^] ^e^l^-'^Z^. (16) 

Since Z = Y.Li N{Li)e~^'^"''"'~^^\ where N{Lj) is the number of dislocation configura- 
tions of total length Lj, the average total dislocation length in the ensemble is 



~ ELN{L)e-^-^ ^^^^ 

i.e., the average total dislocation length is equal to the average number of loops times the 
average loop length. 

The dislocation density, p, is the average total length per unit volume. It then follows 
from (12), (13) and (17) that 

bMT) -^b^-^W/^ e^/'^^^m- (18) 
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3 New melting relation 



The effective line tension, aes, [see Eq. (1)] vanishes at the critical temperature kBTcr = 
ab/ In z'. Consequently, dislocations proliferate as Tcr is approached from below. At tem- 
peratures above Tcr, the divergence of Zi signals the breakdown of the underlying theory, 
and the system enters a new phase. Hence, the temperature Tcr corresponds to a phase 
transition, in which dislocations are copiously produced in the solid. We therefore equate 
the melting temperature, T^, to Tcr- 

The line tension, i.e., the dislocation self-energy per unit length, is assumed to be 
that of a dislocation in a complex array, or tangle, of other dislocations. In that case the 
stress field of a given dislocation beyond ~ -R/2, where R is the mean interdislocation 
separation, is largely cancelled out by the stress fields of the other dislocations in the 
complex array |2^, The line tension is then the sum of the core energy plus the 
elastic energy inside a cylinder of radius R/2 p6| : 



4:71 Kb 2 J Svr V^^V, 

Here, k is 1 for a screw dislocation and (1 — z/)^^ ^ 3/2 for an edge dislocation, i/ being 
the Poisson ratio. Also, G is the shear modulus, b is the Burgers vector magnitude, and 
a is a constant of order unity. In the second half of this equation we have taken distance 
R to be approximately equal to l/^/P; where p is the dislocation density defined in Eq. 
(18). An expression of the form (19) with R = p~^/^ for the dislocation self-energy was 
originally proposed by Mizushima 0|, later put on a sound theoretical basis by Yamamoto 
and Izuyama |l^, and was recently employed by Kierfeld and Vinokur [l^ to model 



dislocation-mediated phase transitions of vortex-line lattices in high-T^ superconductors. 

The constant a accounts for the nonlinear elastic effects in the dislocation core. Hirth 
and Lothe compare dislocation energies in the Peierls-Nabarro (discrete) and Volterra 
(continuum) dislocation models and find 



1 d j sin P 



a eb VeTfl 



+ cos' (3 , (20) 



where 7 = (1 — 2z/)/4(l — i/) ^ 1/8, d is the interplanar spacing, and /3 is the angle between 
the Burgers and sense vectors of the dislocation. In a face-centered cubic (fee) crystal, 
the smallest perfect-dislocation Burgers vectors are i(110)a, and the primary glide planes 
are {111} with d = a/\/3, where a is the lattice constant. Experimental evidence (ref. 
p6| . Table 9-2, p. 275) suggests that the predominant high-temperature glide system in 



body-centered cubic (bcc) lattices is {110}, which has d = a/\/2. The smallest bcc perfect- 
dislocation Burgers vectors are |(lll)a. Thus, in both cases d/b = \/2/3. Averaging over 



P, we find a ~ 2.9 for both fee and bcc lattices. Atomistic calculations of core energies in 
ionic crystals (ref. ||2^, p. 232) indicate that a ~ 3. In metals, no such calculations have 
been performed. We use a = 2.9 for all elements. 

We have also assumed that no backtracking is allowed for dislocations, z' = z — 1, 
since each backtracking would result in a divergence in the linear elastic interaction energy 
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between the overlapping segments. The coordination numbers for the elements considered 
in our analysis below are z = Q for a simple cubic (sc) lattice, z = 8 for bcc and body- 
centered tetragonal (bet) lattices, and z = 12 for fee, hexagonal close-packed (hep), and 
double hep (dhcp) lattices. Replacing 

= Xvws, 

where vws is the volume of the Wigner-Seitz cell of the crystal lattice and A is a geometric 
constant, we finally obtain our formula for the melting temperature of the elements: 

XGvws r-i _ , / 1-45 
Tm= , . , , 7T, S =K\n\^^^=\. (21) 



In ref. ||T9|] we evaluated Gvws/'^'^Tm\n{z — 1) for 51 elements and found 6/X to be 
1.01±0.17, where the error is the root-mean-square deviation. These 6/\ are summarized 
in Fig. 1. 
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Fig. 1. Values of S/X = Gvws /4:'n'Tmln{z — 1) from experimental data for 51 elements. 
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We now assume that the dislocation ensemble is dominated by perfect dislocations 
with the smallest possible Burgers vectors, since dislocation energy is proportional to 6^. 
For a bcc crystal b = av^/2, vws = and b = a/\/2, vws = for a foe crystal. 

Then, b^ ^ 1.30 vws and lAlvws, respectively. For a hep lattice, A = (4/v^) (c/a)~^. 



so that for an ideal hep crystal (c/a = yS/S) one would have A = y/2. As estimated 
for two hep metals [|T3|, b^ ~ 1.42 vws for Mg and 1.24 vws for Zn. Hence, we take 
A = 1.33 ± 0.09 ~ 4/3. This embraces all of the values quoted above. 

In an ensemble of loops there are roughly equal amounts of edge and screw dislocation 
in the crystal, so we have 1/k = (1 — i//2) ± z//2 ~ 5/6 ± 1/6. Therefore, as follows from 
(21), 

/ 2.1 \ 2 (5/6 ±1/6) 



In 



\b'^p{Tm)J (1.33 ±0.09)(1.01± 0.17) 



1.24 ±0.33. 



(22) 



Hence, 



p(T^) = (0.61 ±0.20) (23) 
It follows from Eqs. (21)-(23), with kX = 1.6 ± 0.3, that to ~ 20% accuracy 

Gvws 



^"^ 47rln(;2-l)' 
We regard Eq. (24) as a new dislocation melting law. 



(24) 



4 Free energy of the dislocation ensemble 

To calculate the free energy of a dislocation ensemble, Eq. (15), let us rewrite Eqs. (10) 
and (12), using Eq. (19) with R = Ij and replace the sums (which start with L = 46, 
the smallest loop length) by the corresponding integrals: 



Zi = -ri^A(q, z ) / dxx~ 

0"* J4 



/■oo [Ah'^n 
i{T) = A{q,z') / dxx-'^ -^{^'f 



K 



Gb^ 



(25) 



(26) 



Here, {Ab'^p/a'^) {z'y/" = exp{— 87ro-cfr/'^G'6^} < 1, since a^s > 0. Integrating Eq. (25) by 
parts we find 



V I KGb' , 

_ V_ 
qb^ 



A{q,z') 
41 



4c N 



Hence, 



F = -kBTe^^'^^Z^ = ^ (^aespb^ - 



A{q,z')e 



'\pfJ./kBT 



41 



(27) 



:28) 
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where we have replaced (nGb^/STr) ln(a;^/46^p) — fcsTln^;' by bacs, in view of (1),(19), 
and used Eq. (18). 

The second term on the right-hand side of Eq. (28) takes its largest value at T = T^, 
where cTefr = 0. To estimate its contribution to the free energy, consider the case of Cu 
discussed in more detail below. In this case, to estimate A{q, z') exp{fi(Tm)/kBTm}/4:'^, 
we use Eqs. (12), (18), and replace the sum by an integral: 

A(q, z')e^^(T^y'^T^ _ A{q, z') b'pjTj _ b'pjTj _ b'pjTj {q - 1) 
4-? 4-? ^{Tm) 4iJ^dxx-'i 4 

As discussed in Section 2, the value of q may be expected between 3/2 (Brownian loops) 
and ~ 7/4 (self-avoiding loops). With b'^piTm) given in Eq. (23), we therefore obtain 
A(g, z') exp{/i(T^)/fcBT^}/4« = 0.095 ± 0.037 ^ 0.1. 

o 

Hence, the contribution of the second term to qF/V would be ~ —0.7 meV A^^. As 
seen in Fig. 2, this contribution is negligibly small. In fact, the second zero of F for 
T = Tm and 74(g, z') exp{/i(Tm)/fc_BTm}/4'' = 0.1 occurs at fe^p = 0.61, which is within 
~ 5% of the value of 0.64 (the second zero of F at T = with F/V given in (29)), and 
within uncertainties in the values of b'^p{Tm) in Eq. (23). 

Thus, we have derived the dislocation free energy density, and it is given approximately 

by 

^ ^ ~ f^effP = - — ^ +^ln(z-l) UV (29) 



V ^ \M \a'' ) b^ 

This form for the free energy density was previously suggested but not derived by Cotterill 



15 1 . It was later put on a firm theoretical basis by Yamamoto and Izuyama [|T^. It also is 
a fundamental ingredient in the recently developed theory of dislocation-mediated phase 
transitions of vortex-line lattices in high-T^ superconductors |18| . 



In Fig. 2 we plot qF(p)/V from (29) for Cu for three different temperatures: T <T, 



mi 



T = Tm andT > T^. We take k = 6/5, a = 2.9, G = 47.7 GPa |8|, = 1356 

o 

K, and b = 2.55 A [|19|. A first-order phase transition, that is melting, takes place 
when the second zero of F{p) occurs at the critical dislocation density, p(Tm). This is a 
transition from a perfect crystalline solid to a highly dislocated solid, not a liquid. In 
fact, our theory describes dislocations, which do not exist in liquids. If a dislocation is 
viewed as a disclination dipole p9|, the dislocated solid may in turn undergo a Kosterlitz- 



Thouless-like transition [g^l to a phase of free disclinations, i.e., a liquid. This dislocated 
solid may then be viewed as the three-dimensional analog of an intermediate hexatic 
phase, between a solid and a liquid, in the Halperin-Nelson theory of two-dimensional 



melting |31|]. The clarification of this point needs further investigation, to be undertaken 



elsewhere. Patashinskii et al. |]TT| also identified melting as a transition from a perfect 



crystalline solid to a highly dislocated solid, and Nelson and Toner found residual 



bond-orientational order in a three-dimensional solid with an equilibrium concentration 
of unbound dislocation loops, which is analogous to that in the two-dimensional hexatic 
phase. 

Note that it is not possible to increase the dislocation density progressively from zero 
to p{Tjn) at a temperature lower than (e.g., by deformation) because of the high 
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energy barrier at the maximum of qF{p)/V. Hence, the dislocation density, as a function 
of temperature, is 



PiT) 



0, 



T < T„ 



p{Tm), T = Tm- 



(30) 



In fact, it can be shown that (30) is the only physical solution of (12) written as a "gap'' 
equation: 



n=4 



KGPn/SnkBT 




5 Latent heat of fusion 

For the ensemble of strings on a lattice considered in Section 2, the internal energy and 
pressure are 
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P = kBT 
Hence, the enthalpy is 



dV 



H = U + PV 



y_ 

63' 



,illkBT 



ahiiT) + ksT liT) 



(32) 



(33) 



The latent heat of fusion is the enthalpy difference: 

L^ = H{TJ-H{0). (34) 

In our case, H{0) = 0, which follows directly from (30)-(33) and Eq. (18). Using <Jcs(Tm) = 
and the melting condition ksTm = crb/ \n{z — 1), we obtain 

Lm = ^e^^'^-^/'^^^-kBT^ [Hz - l)e(r^) + e(T™)] . (35) 

To obtain the latent heat per mole, the quantity tabulated in the literature, one has to 
multiply the expression (35) by the ratio of the number of atoms per mole, A^^, to the 
total number of atoms in the volume V, which is equal to V/vws- Replacing N^kB by the 
gas constant R, and using = Xvws, we obtain 



A 



■e(Tji?T„ln(z-l) 



1 liT„ 

\n{z - 1) e(r„ 



(36) 



To estimate the ratio ^(Tm)/^(Tm), we replace the sums in Eqs. (12), (13) by the 
corresponding integrals: 



e(T^) _ /la=4 d{Llh) (L/b) 



-1-1 



q-1 



aTm) jr/b=Ad{L/b){L/b)-<^ Aq ' 
With 3/2 < g < 7/4, as discussed in Section 2, 0.083 < (g - l)/4g < 0.107, i.e., 

e(T™ 



(37) 



0.095 ±0.012. 



(3^ 



Therefore, the contribution of the second bracketed term on the right-hand side of Eq. 
(36) (corresponding to the work contribution to enthalpy) is 0.04 — 0.06 {Q < z < 12). 
We expect, therefore, that neglecting the second bracketed term on the right-hand side 
of Eq. (36) will introduce an error not larger than ~ 6%. Hence, with accuracy of ~ 94% 
we have the following formula for the latent heats of the elements: 



A 



b'^p{Tm)RTm ln{z 



(39) 



where we have replaced e^^'^"'^^''^'^"'^{Tm) by b'^p{Tm), in view of (18). The proportionality 
of latent heat of fusion to the critical concentration of defects (multiplied by the core 



energy) has been noted previously by Cotterill ||T5 
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In Fig. 3 we plot the values of b'^p{Tm) extracted from the experimental data on latent 
heats for 75 elements. For this analysis, the values of both Tm and are mostly taken 
from ||33|. For Be, Hf, Sc, Sr, Y, the lanthanides Dy, Ce, Er, Gd, Ho, La, Nd, Sm, Tb, Yb, 



and the actinides Am, Cm, Th, we disregard their high-T bcc phases which exist only in 
the very vicinity of melting. (The intermediate hep— s>fcc phase transition for Yb, dhcp— 
fee for Am, Ce and La, and fee— >hcp for Sr, as well as hcp^fcc for Co, do not change 
coordination number.) The crystal structure chosen for the evaluation of Ca, Co, Mn, N, 
Np, O, Sm, Ti, Tl, U and Zr corresponds to the phase from which melting occurs. The 
data on both T, 
and L 



and Lm for H, N, O, Pa and Rn are taken from |34 
for Am and Cm, and on L 



for Ar, Kr, Ne and Xe are taken from 35 



The data on both 
The 



T 

data on Lm for the lanthanides are taken from The following values of A are used: 1 
for sc, 1.3 for bcc, 1.41 for fee, 1.24 for Zn, 1.42 for Mg, and 1.33 for all other elements. 
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Fig. 3. Critical dislocation density as extracted from the experimental data on latent heat 
of fusion for 75 elements. 
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For all these elements we find 



p(T^) = (0.66 ±0.11) (40) 

where the error is the root-mean-square deviation. This value is in good agreement with 
that obtained from the melting temperatures alone, Eq. (23). 

Note that the possible inaccuracy in the value of A for the hep, dhcp and bet elements 
used in this analysis, on the order of ~ 7%, may slightly increase uncertainty in the value 
of p(T„) in Eqs. (23) and (40). 

We do not have a reasonable explanation for the anomalously high values of p{Tm) 
for the noble gases. If the noble gases are excluded from the analysis, then p{Tm) turns 
out to be (0.63 ± 0.06) b'"^ for the remaining 70 elements. Note also that for deuterium 



(D), which is not included in Fig. 3, with the data on and Lm from ||3^, we obtain 
b^piTm) = 0.70. 

The uncertainty-weighted average of the values of p{Tm) given in Eqs. (23) and (40) 



IS 



p(Tj = (0.64 ±0.14) 6-2, (41) 
which we take as our result for the critical dislocation density at melt. 



6 Volume change at melt 



As an independent consistency check on the relations (23) and (40), we determine the 



critical dislocation density using the formula |13, 37| 



V 



2^5 I 



77 b'piTn 



(42) 



where AV^ is the difference between the liquid and solid specific volumes at melt, G and 
B are the shear and bulk moduli, respectively, and 70 is the Griineisen constant. Here, 
e is identified with the dilation of the lattice as the reaction of the crystal to the sudden 
proliferation of dislocations. In Table 1 we show the values of h'^p{Tm) calculated for 32 
elements for which we could find zero-pressure data on 7^ and e. The experimental values 
of e are mostly taken from ref. pSl, and those of G, B and 7g from For Ar, Kr, Ne 



and Xe, the values of G and B are taken from [^, and those of 7^ from 
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element 


B, GPa 


G, GPa 


7G 


e 


h^p{Tm) from Eq. (42) 


h-^p{Tm) from Eq. (39) 


Ag 


103 


29.8 


2.40 


0.052 


0.39 


0.67 


Al 


76.0 


26.1 


2.19 


0.064 


0.44 


0.81 


Ar 


1.83 


0.75 


2.59 


0.144 


0.69 


1.00 


Au 


173 


28.0 


2.99 


0.055 


0.57 


0.65 


Be 


111 


151 


1.11 


0.115 


0.51 


0.63 


Ca 


16.7 


7.4 


1.15 


0.048 


0.64 


0.63 


Cs 


2.01 


0.65 


1.41 


0.026 


0.36 


0.56 


Cu 


137 


47.7 


2.02 


0.046 


0.35 


0.68 


Eu 


17.0 


7.53 


1.39 


0.048 


0.50 


0.67 


Gd 


37.8 


21.6 


0.63 


0.021 


0.59 


0.65 


Ho 


40.8 


26.3 


1.18 


0.075 


0.65 


0.66 


In 


42.0 


4.78 


2.43 


0.025 


0.49 


0.62 


K 


3.3 


0.9 


1.29 


0.025 


0.46 


0.56 


Kr 


2.04 


0.85 


2.64 


0.151 


0.70 


1.00 


Li 


12.1 


3.85 


0.92 


0.016 


0.41 


0.53 


Lu 


47.6 


27.2 


1.06 


0.036 


0.41 


0.67 


Na 


6.74 


1.98 


1.19 


0.027 


0.52 


0.56 


Nb 


171 


37.6 


1.77 


0.029 


0.44 


0.59 


Nd 


32.9 


17.4 


0.57 


0.009 


0.34 


0.56 


Nc 


0.88 


0.40 


2.79 


0.156 


0.62 


0.98 


Ni 


183 


85.8 


1.93 


0.063 


0.37 


0.72 


Pb 


44.7 


8.6 


2.74 


0.037 


0.36 


0.56 


Pd 


193 


48.0 


2.56 


0.059 


0.47 


0.66 


Pt 


283 


63.7 


2.87 


0.066 


0.51 


0.68 


Rb 


2.3 


0.63 


0.99 


0.026 


0.70 


0.60 


Ta 


193 


69.0 


1.74 


0.052 


0.50 


0.59 


ib 


38.7 


22.1 


0.74 


0.032 


0.65 


0.67 


Tl 


35.7 


5.4 


2.10 


0.033 


0.60 


0.60 


Tni 


46.2 


29.1 


1.43 


0.069 


0.47 


0.68 


W 


310 


160 


1.67 


0.090 


0.63 


0.77 


Xe 


2.1 


1.0 


2.56 


0.130 


0.54 


1.01 


Yb 


14.9 


8.06 


1.04 


0.036 


0.44 


0.56 



Table 1. Values of lP'p{Tm) from experimental data on volume change at melt for 32 ele- 
ments. For comparison, we also show values of h^p{Tm) extracted for the same elements 
from the data on latent heats. 



For all 32 elements in Table 1 we find 

= (0.51 ± 0.11) r^ (43) 

where the error is the root-mean-square deviation. This is somewhat lower than but still 
in agreement with both Eqs. (23) and (40) taking into account uncertainties associated 
with the three values. 

For comparison, we show in the last column of Table 1 the values of b^p{Tm) extracted 
for the same elements from the data on L^. It is seen that the agreement between two 
sets of the values of 6^p(T^) is reasonably good, except for Ag, Al, Cs, Cu, Ni, Pb and 
Xe, for which the difference in both values of b'^p{Tm) is on the order of ~ 50 — 60%, Lu, 
Ne and Nd for which the difference is ~ 45%, and Ar, Kr, Pd and Tm, for which it is 
~ 35%. For all other elements, the difference does not exceed ~ 30%. 
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Note that the contribution of the volume change at melt, e, to the latent heat of 
fusion is proportional to -C 1 [|1^, 0, and is therefore negligibly small compared to 
the right-hand side of Eq. (39). 



7 Concluding remarks 

Our theory of dislocation-mediated melting was developed in the approximation that 
dislocations are non-interacting. This approximation is good only in the vicinity of melt 
where the dislocation density is very high and the otherwise long-range interactions are 
sufficiently screened [^. The statistical mechanics of non- interacting dislocations on 



a lattice yields simple, accurate relations between the dislocation density at melt and 
both the melting temperature and latent heat of fusion, despite the indeterminacy of the 
parameter q that takes into account the possible non-Brownian nature of the dislocation 
network. The values of p{Tm), as determined from an extensive analysis of Tm and Lm 
data, are remarkably consistent: (0.61 ± 0.20) and (0.66 ± 0.11) respectively. 
The uncertainty-weighted average of these values is p{Tm) = (0.64 ± 0.14)6"^, which we 
take as our result for the dislocation density at melt. Poirier and Price [|13| analyzed 
14 elements and found p{Tm)vws/b = 0.48 ± 0.12. Using vws = b^/^ with A ^ 4/3, 
their result corresponds to p{Tm) = (0.64 ± 0.16) which is in excellent agreement 
with ours. Kierfeld and Vinokur |jl8| modelled dislocation-mediated phase transitions of 
a vortex-line lattice and found piT^) ^ 0.6 Vachaspati's |HOi study of topological 



defect formation gave a?p{Tm) ~ 0.88 for a simple cubic lattice. This translates into 
p(T„) ^ 0.66 for bcc lattices (a = 2/^3 h) and p(r„) ^ 0.44 fo-^ for fee lattices 
(a = a/2&), which are consistent with our result. In agreement with Vachaspati, Kibble 
PT| found a^piTm) ~ 0.89 for a simple cubic lattice. 

Although our main results do not depend on the precise value of g, there is a particular 
value of q at which the relations (29) and (39) become exact: g = 1. In this limit, as seen 
in (12), iiTm) — ^ oo, so that Eq. (39) becomes exact in view of (36). Requiring finite 
internal energy in this limit leads, via (31), to exp{/i(Tm)//cB7'm} — > {piT^) —oo), 
and therefore, Eq. (29) becomes exact, since the second term on the right-hand side of (28) 
disappears. In fact, the study of cosmological networks of string loops in 3 dimensions 
by Magueijo, Sandvik and Steer [H2| results in a scale-invariant loop distribution of the 



form of Eqs. (9), (10) with g + 1 < 5/2 : [|| 1.9 < g + 1 < 2.1, or [|| g + 1 = 2.03 



(plus error bars), and so in this study g ~ 1. Thus, it is quite possible that linear defects 
which correspond to two apparently distinct physical phenomena, namely cosmic strings 
and crystal dislocations, are of a very similar statistical-mechanical nature. 

The average total dislocation length per Wigner-Seitz cell at melt is piTm)vws = 
bp{Tm)/\ ~ b/2, since A ~ 4/3. Since a Wigner-Seitz cell contains z links, each of length 
6/2, it follows that, on average, one of z links in each Wigner-Seitz cell is covered by a 
dislocation. Since each such a link is shared between two atoms, on average, half of the 
atoms are within a dislocation core at melt. 

If we use p{Tm) = 0.64 fe-^, then to ~ 20% accuracy the melting temperatures and 
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latent heats are given by 

''^- = l^Y '''' 

L„ = iH(i_li)iiT„. (45) 

The accuracy of these relations depends critically on the factor of ln(2; — 1), which is 
characteristic of a theory based on line-like degrees of freedom. 
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